DFPD 02/TH/03 
hep-th/ 02020821 

February 2002 



The beta function of N=l SYM 
in Differential Renormalization 



Javier MAS a 'Q, Manuel Perez- Victoria b f] and Cesar SEUAS a £] 

a Departamento de Fisica de Particulas, Universidade de Santiago de Compostela, 

E- 15706 Santiago de Compostela, Spain 

b Dipartimento di Fisica "G. Galilei" , Universitd di Padova and 
INFN, Sezione di Padova, Via Marzolo 8, 1-35131 Padua, Italy 



Abstract 

Using differential renormalization, we calculate the complete two-point function of the 
background gauge superfield in pure N=l supersymmetric Yang-Mills theory to two loops. 
Ultraviolet and (off-shell) infrared divergences are renormalized in position and momen- 
tum space respectively. This allows us to reobtain the beta function from the dependence 
on the ultraviolet renormalization scale in an infrared-safe way. The two-loop coefficient 
of the beta function is generated by the one-loop ultraviolet renormalization of the quan- 
tum gauge field via nonlocal terms which are infrared divergent on shell. We also discuss 
the connection of the beta function to the flow of the Wilsonian coupling. 
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1 Introduction 



In this paper we revisit a somewhat old controversy: the origin of higher-order perturba- 
tive contributions to the beta function in supersymmetric gauge theories. The relevance 
of infrared (IR) modes and the distinction between the Wilsonian action and the gen- 
erator of 1PI functions are at the heart of this debate. A good understanding of these 
issues is relevant, for instance, for the comparison of field theory results with holographic 
renormalization group flows 

The so-called "exact beta function" of general N=l SYM was discovered by Novikov, 
Shifman, Vainshtein and Zakharov (NSVZ) using instanton analysis |§. In this approach 
it is clear that corrections to the one-loop result have an IR origin in an imbalance in 
the number of fermionic and bosonic zero modes. For pure SYM the NSVZ beta function 
reads 

(3(9) NSV z = ig 2 ] cTO* ■ 

This formula was first derived by Jones in ||. Later on, superspace perturbative com- 
putations to two loops were carried out using regularization by dimensional reduction 
(DRed) 0, |5|, p. In this method the two-loop correction to the beta function arises from 
a local evanescent operator specific to DRed|. This operator is not available in regulariza- 
tion schemes that stay in four dimensions and Grisaru, Milewski and Zanon pointed out 
that this seems to imply that no divergence should occur beyond one loop, in conflict with 
the DRed result || . NSVZ then observed that in four dimensions the higher-loop correc- 
tions can only arise from nonlocal operators that are nonanalytic at vanishing external 
momentum j7|. This behaviour can appear only from the domain of virtual momenta of 
order of the external momenta. Since this domain is excluded by definition in the Wilson 
action, the flow of the Wilsonian coupling constant is purely one-loop ||. The standard 
modern proof of this fact is based on the holomorphic dependence on the complexified 
coupling constant. The running of the physical coupling constant, on the other hand, 
has higher-order contributions that appear when one takes the expectation value (in an 

1 The calculation of the beta function of N=l SYM in the context of the AdS/CFT correspondence 
has been recently considered in jl| . In order to discriminate between the different possible results (which 
depend on the detailed UV/IR relation one uses), a prcliminar requisite is to identify the field theory 
answer one would like to reproduce. This requires a correct field-theoretical interpretation of the role of 
the IR degrees of freedom. 

2 This is the only local gauge invariant operator of the appropriate dimension: Tr J d i xd i 6T a Tb(8 b a — 8 b a ), 
where S and 5 are Kronecker deltas in four and n dimensions, respectively. Using a Bianchi identity it 
can be cast in a form proportional to the classical action: — eTr J d 4 xd 2 8W a W a - 
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external field) of the operators in the Wilson action. Furthermore, the IR pole is related 
to an anomaly, and this is the crucial fact that allows to determine exactly the higher 
order contributions [[J. In a later contribution, Arkani-Hamed and Murayama rederived 
the NS VZ beta function in a purely Wilsonian setup . They showed that the canonical 
Wilsonian coupling constant obeys a NSVZ flow. The reason is that keeping canonical 
kinetic terms at each scale requires a rescaling of gauge field which is anomalous. This 
anomaly generates the higher order corrections. These authors claimed that their calcu- 
lation only depends on ultraviolet (UV) properties of the theory, and thus questioned the 
IR origin of the corrections. Furthermore, they pointed out that the method of differential 
renormalization |E2|] clearly displays an UV origin of the corrections. This interpretation 
has been criticized in |10f . 

The purpose of this paper is to try to contribute to the understanding of these issues by 
performing an explicit calculation in differential renormalization (DiffR) ||12|| . The interest 
of using this method is twofold: on the one hand we are able to derive the beta function di- 
rectly from the scale dependence of finite renormalized Green functions, rather than from 
"infinite" counterterms; in this derivation we see explicitly how nonlocal terms contribute 
to the beta function in a perturbative four- dimensional method. On the other hand, we 
clearly separate UV divergences from the off-shell IR divergences that afflict these calcu- 
lationsQ. As a by-product, we develop some calculational tools that we believe have an 
intrinsic interest to the SUSY-community. In fact, DiffR is a computational program that 
seems to be especially taylored for supersymmetric theories: it neither requires continua- 
tion in spacetime dimensions, nor changing the field content of the bare lagrangian. It is 
rather an implementation of Bogoliubov's R operation, which yields directly renormalized 
correlation functions satisfying renormalization group equations. DiffR has been applied 
before to supersymmetric calculations in the Wess-Zumino model Jl3| to three loops, in 
pure SYM and SQED to one and two loops, respectively |14|], and in supergravity to 
one loop |T5|. Its implementation in symbolic programs (TjJ enables efficient one-loop 
calculations in more involved models like the MSSM |T7 . 



The layout of the paper is the following. In Section ||] we review the method of 
DiffR and introduce the new tools needed for our calculation. In Section |3| we quickly 
review the supercovariant background field method to settle the notation. In Section |] 
we calculate, to two loops, the renormalized correlation function of two background gauge 

3 In the dimensional methods UV and IR divergences are mixed. In H) the IR divergences were removed 
by the choice of a nonlocal gauge fixing that kept the renormalized quantum propagator in the Feynman 
gauge. In g the IR divergences were simply subtracted by an R operation pl[ . We shall also perfom an 
R operation, but keeping track of the resulting finite part. 
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superfields in supersymmetric gluodynamics. The contribution of matter fields can also be 
computed using the techniques described here, but the results will be presented elsewhere. 
In Section |5] we write the renormalization group equation and determine the first two 
coefficients of the beta function. Section ^ is devoted to a discussion of the origin of the 
higher-order coefficients in our calculation and in previous works. In particular, we discuss 
the relation between our result and the flow of the Wilsonian coupling constant. Finally, 
in the Appendix we compute the two-point function of quantum gauge superfields at one 
loop, and determine the renormalization group coefficient associated to the gauge-fixing 
parameter. 



2 Differential renormalization 



Differential renormalization |T2[ is a method that defines renormalized correlation func- 
tions without intermediate regulator or counterterms. This is achieved by writing coor- 
dinate-space expressions that are too singular as derivatives of less singular ones. The 
derivatives are then understood in the sense of distribution theory, i.e., they are prescribed 
to act formally by parts on test functions, neglecting divergent surface terms. Diagrams 
with subdivergences are renormalized according to Bogoliubov's recursion formula. This 
procedure leads to bona fide distributions that respect the requirements of quantum field 
theory. Consider as a simple example the singular function (1/x 2 ) 2 . Its renormalized form 
is simply expressed by 
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1 lnx 2 M 2 . . . 

= -tO 5 — + a uv 5(x), 2.1 

R 4 X 



where the D'Alembertian acts "by parts" . Note that the bare and the renormalized expres- 
sions, when understood as functions, coincide for i^O. However, only the renormalized 
expression is a finite distribution on test functions defined over the complete space. The 
arbitrary scale M, with dimension of mass, must be included for dimensional reasons and 
plays the role of a renormalization group scale. Related to this is the fact that one can al- 
ways add a local term of the appropriate dimension, which reflects the freedom in choosing 
a scheme. Note that the arbitrary constant ajjy can be absorbed into a redefinition of M. 
Different renormalized expressions can in principle have different renormalization scales 
and/or different local terms. Our approach will be to write a single (UV) renormalization 
scale and adjust the contact terms in such a way that gauge invariance is preserved. 
Analogously, IR divergent expressions can be made finite by differential renormaliza- 
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tion in momentum space. For instance, 



1 

P 1 



R 



\np 2 /M 2 IR 



P 



+ a IR 5(p) . 



(2.2) 



We have denned for convenience Mm = IMmllE-, where 7^ is Euler's constant, and 
distinguished the IR scale from the UV one. This is an explicit realization of the so-called 
R operation that subtracts IR divergences. Again, diagrams with IR sub divergences are 



treated according to a recursion formula |TT| , |18fl analogous to the UV one. 

Since UV and IR overall divergences are local in coordinate and momentum space, 
respectively, the R and R operations commute, and one can define an operation R* = RR 
to renormalize both UV and IR divergences [O, [TBI]. The fact that the UV and IR 



renormalizations decouple means that the UV and IR renormalization scales should be 
independent. In DiflR this can be achieved by a careful adjustment of the local terms 
involving both scales^. Let us implement this idea in an example which will play a central 
role in the calculation of Section HI. Consider the IR singular expression 



In p 2 /M 2 
p 1 



(2.3) 



that arises in IR divergent expressions after renormalization of a UV subdivergence (with 
ajjv=0)- The consistent IR renormalization of (|2.3|) is given by 



\np 2 /M 2 
p~ 4 



1 



A" 



R 



\n 2 p 2 /M 2 R + 2 \np 2 /M 2 R (1 + \np 2 /M' 

p2 



+ (a IR In 



— ^ + b 

M 2 



(2.4) 



This expression differs from the usual one given in by scale- dependent local terms 
proportional to In 2 M 2 /Mf R (appart from the explicit local terms with coefficients am and 
b IR ). It should be used whenever the "new" scale is to be treated as independent from the 
"old" one, for consistency of the loop expansion. The scale-dependent local terms of ( |2.4| ) 
are fixed by the requirement that the IR renormalization commute with a rescaling of M, 
that is to say, 



M 



d 
dM 



\np 2 /M 2 
p 1 



M 



d \np 2 /M 2 



dM 



P 



(2.5) 



R 



4 IR DiffR was investigated in [[19| where it was concluded that the combination of UV and IR DiflfR 
was inconsistent, as the results depended on the order in which integrations were performed. According 
to p0| , however, this corresponds to the natural arbitrariness of the IR renormalization, and this author 
has actually proposed in a consistent version of DiffR that deals with both UV and IR divergences. 
Our approach here will be closer to the original version of DiffR 
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Observe that the UV scale M only appears in ( f2.4| ) in single logarithms. This is fine, for 
double logarithms of M are expected to appear only when the bare expression contains 
both a UV subdivergence and a UV overall divergence. Observe also that a rescaling 
of Mir in (|2.4J) gives a local term in p-space. This procedure can be extended to more 
general situations, but the identity Q2.4Q is all we need for the calculation at hand. 

Let us finally deal with the scale-independent local terms. They must be chosen such 
that the renormalized correlation functions respect the fundamental symmetries of the 
theory. In our problem these are supersymmetry and gauge invariance. Since the first 
is automatically preserved in superspace, we only have to worry about the second. In 
ordinary DiffR, we would have to study the Ward identities order by order and adjust the 
local terms by hand so that they are satisfied. For the calculation of the two-loop beta 
function it is sufficient to impose the Ward identities at the one-loop level, but this is 
complicated in the framework of covariant supergraphs. Life gets much easier, however, 
when one uses the so-called constrained differential renormalization (CDR) ||23|| . This is 
a procedure that fixes the arbitrary local terms a priori in such a way that the Ward 
identities are directly fulfilled |23|, |24], p5|] . Furthermore, CDR respects supersymmetry 
in component field calculations |L5], |2~6fl . In superspace calculations CDR is particularly 
simple because, after performing the superalgebra, all sub divergences are Lorentz scalars. 
According to the CDR prescriptions, this means that the local terms in the renormalized 
subdiagrams are universal, i.e., they are independent of the Green function or diagram 
they appear in. Specifically, in our calculation we shall take cljjv = am = 0, and keep bm 
arbitrary (but unique) till the very end. When calculating beta functions, this is all we 
need. Nevertheless, in order to calculate the complete two-point function, we shall also fix 
by hand the local terms that appear in superficially divergent tensor structures, so that 
gauge invariance be preserved. (We need this straightforward adjustment because CDR 
has not been developed beyond the one loop level.) 



3 Super co variant background field method 

In the background field method, the total gauge superfield Vr is splitted into background 
B and quantum V superfields according to 

e^ = eW ; e B = e"e" (3.1) 
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We will follow closely the notation and conventions given in [22]. Background covariant 
derivatives can be defined as follows 



V a = e~ n D a e n = D a - iT a (tl) 

Va = e A Aie- A = - it A (0) (3.2) 

In the chiral representation for covariant derivatives, the "quantum-background" splitting 
amounts to 

V Q = e- VT D a e VT = e- A e- 9V V a e 9V 

V« = = e _A V d e n . (3.3) 
The classical action of pure N=l SYM then adopts the following form 

Sym = j 2 J d 4 xd 2 6TrW 2 = 7^ J d 4 xd 2 6Tr [V d , {V A , V Q }]) 

= -^|rf 4 x^TrQ[v Q {V Q ,e-^V Q e^}]y (3.4) 

The quantum-gauge fixing V 2 V = retains background covariance. Usual averaging 
requires the introduction of Nielsen-Kallosh ghosts: 

J VfVfVbVb 5(V 2 V - f) 5(V 2 V - /) exp d 4 xd 4 9(ff + bb)^j 

= J VbVb exp(^-^ J d 4 xd 4 9[W(V 2 V 2 + V 2 V 2 )V + bb]^j (3.5) 

Expanding Sym + S g f in powers of the quantum field V yields 

S = Sq + S2 + Dint 

= \J d 4 xd 2 6 TrW 2 - d 4 xd 2 6V[n + £(V 2 V 2 + V 2 V 2 )]V 



- J d 4 xd 4 6 (c'c + cc' + (1 + + S int (V,c,d) 



(3.6) 



where £ = J-l, □ = n-iT a "d a a-{{d a j: a& )-\T a *T a «-W a V a -W«V « and the dots 
stand for terms with higher powers of V or c, d . All anticommuting superfields, c, d and 
b interact with the background field T through the constraint that they be background 
covariantly chiral, Vc = V6 = Vc = V5 = 0. The Effective Action in the Background 
Field Gauge admits a gauge invariant expansion in the form 

T[B] = J d 4 xd 4 yd 2 9 [W a (x)W a (y)] G {2 \y - x) + ... (3.7) 
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Our aim is to calculate the 2-point 1PI function G^ 2 \ For the perturbative computation 
of T, we expand the action in powers of V and distinguish the "free" part (in the presence 
of the background), So, from the interacting part, Si nt : 

T[B] = s + rl loop 

"AAA 

5 J' 5j' 5j 



+ exp S- n 



exp 



dvei-Jn J-jn+j 



(3- 



J=j=j=0 

|(V a W af )]$ is the D'Alembertian acting on 



where = V 2 V 2 $ = [□ - iW a V a 
background-covariant chiral fields (similarly for anti-chiral fields = V 2 V 2 $ = [□ — 
zW"V A -f (V Wd)]$), and J, j,j collectivelly denote sources for vector chiral superfields 
in ( |3.6| ). r^ oop stands for the 1-loop contribution in the gauge £. We are interested in 
computing the two-point amplitude at two loops in the Feynman gauge £ = 0. However 
the gauge parameter £ is renormalized and the RG equation will generically contain a 
term ^^d/d^G^ |^ =0 . Since 7^ = <9£/<91ogM = 0{g 2 ) + ... the linear dependence in £ of 
the 1-loop Green's function will be needed. 



4 Calculation of the two-point function 



4.1 One Loop 

Formally, the exact one-loop contribution in an arbitrary Lorentz gauge, £, is 



-plloop 



Air In (a + £(V 2 V 2 + V 2 V 2 )) + 2Tr ln( V 2 V 2 ) + Tr In ((1 + f ) V 2 V 2 ) 



1.1) 



The three terms represent a loop of quantum gauge superfields, Faddeev-Popov ghosts 
and Nielsen-Kallosh ghosts, respectively. Expanding to linear order in £ 



illoop 



-iTrlnn + 3Trln(V 2 V 2 ) + |tt 



— V 2 V 2 + h.c. 



+ 0{i 2 ) . (4.2) 



From this expression, we are instructed to expand in powers of the external background 
field ft. The first term in (|4.2|) is well known to start contributing from four-point functions 
up ||22|| . The second piece, stemming from the ghosts, yields the standard contribution 
to the 1-loop beta function in the Feynman gauge. The diagram involved is shown in 



Figure 1-a) and yields JTJ 

3C A 



d A xd 4 yd A 6W a (x)T a (y)A 2 (x - y) 



R 



3C A 
4% 2 



d 4 xd A yd 2 9W a (x)W a (y)n A(x - y) ln(x - y) 2 M 2 (4.3) 
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w / \ r w / ^ r 

y\AA /WV, 7W\ /VAA, 



a) 



b) 



Figure 1: Diagrams contributing to Y 



Hoop 



after use of (|2.1|) with auv = 0. Here the free propagator is A(x) = — (x 



□ 



1 1 

47T 2 X 2 

The last term in (|4.2| ) involves corrections to the gauge parameter. Standard V-algebra 
manipulations reduce it to 



<Fx<Fyd z e W a (x)W a (y)n 



A(x-y) J d A zA(z - y)A(z - x) 



(4.4) 



The integral that remains to be done corresponds to the insertion diagram of figure lb), 
and diverges logarithmically for large \z\. This is the first instance of an off-shell IR 
singularity, that we renormalize along the lines explained in Section § 



/ 



d 4 z 



(x — z) 2 (z — y)' 



f d 4 p \ e ~ ip{x - y) 
J V 



R 



-ir 2 ln(x-y) 2 M 2 R 



ln-rf^ 



(4.5) 



We have used the identity ( |2.2p with am 
contribution to linear order in £: 



0. In the end we obtain the following 



4% 2 



D[A(x - y) \n(x - y) A Mf R 



(4.6) 



and the full one-loop contribution results in 

C A 



-plloop 



4 3 7T 2 



d 4 xd A yd 2 e W a (x)W Q (y)n 3A(x - y) ln((x - y) 2 M 2 ' 
- £A(x - y) ln((x - y) 2 M 2 R )} + 0(£ 2 ) + . . . 



(4.7) 
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4.2 Two Loops 

We will compute the two-loop contribution in the Feynman gauge, £ = 0. From (|3.8|) , the 
higher-loop contributions come from the expansion of a vacuum diagram with propagators 
and vertices made out of background covariant derivatives (Figure 2). In particular there 
is a single such diagram at two loops, as shown in [|]. After performing the V-algebra 




Figure 2: Diagrams contributing to at two loops. 



we find the following nonvanishing contributions (up to a common factor — 3g 2 C\/2): 

D x = tJd 4 xd A yd A 6W«(x,6)W a (y,6)a».[-A(x-y)d,I (x,y)], 

D 2 = iJd 4 xd A yd*9W & {x,9)W a (y,9)<j^ [2A{x - y)d fi I°{x, y) - d^{A(x - y)I°{x, y))] 

D 3 = ^Jd*xd A yd*9r^x,9)r»(y,6)[-d,d u (AI°(x,y)) (4.8) 

+ 40 M ((ft,A) I°(x, y)) - 4((^^A)/°(x, y))] , 
D 4 = l -f d 4 xd 4 yd 4 9T^x,9)T"(y,9) (i^x, y) - I 2 ^(x, y)) , 

D 5 = \ I d 4 xd 4 yd 4 9 r»(x, 9)T»(y, 9) 5^ (nA(x - y))I°(x, y) . 
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Here and in the following, all derivatives act on x unless otherwise stated. In D§ we have 
replaced 8{x) by DA(x — y). The integrals P are defined as follows: 

I°(x,y) = J d 4 ud 4 v A(x-u)A(y-v)A 2 (u-v) 

I^x, y) = J d 4 ud 4 v A(u - v) A(x - v)d^A(x - u) A(y - v)d v v A(y - u) (4.9) 

I 2 v (x, V) = J rf4 wd 4 t> A{u - v) A(x - u)d*A{x - v) A{y - v)d y u A(y - u). 

As shown above, all expressions but -D4 are obtained from a single integral I , which is 
both UV and IR logarithmically divergent. On the other hand, derivatives of I°(x) are 
just UV divergent. The integrals I 1 and I 2 , which appear in D^, are IR safe as well. 

The strategy now is to look for a renormalized expression for the sum Yl^=i To 
begin with we renormalize the integral 1°. We must first cure the UV sub divergence, 
using the DiflR identity (|2J]) with auv = 0: 



d 4 ud 4 v 



(4vr 2 ) 4 J (x- u) 2 (y - v) 2 (u - v) 4 ' 

1 f d'ud'v , 1 ,., 1 J^-^ 



4(47r 2 ) 4 J (x - u) 2 (y - v) 2 (u - vy 

1 f l4 1 ln(a; - v) 2 M 2 
' d v- 



2 



4(47T 2 ) 3 J (y-v) 2 (x-v) 
1 f ^ ln P 2 /M 2 e _ ip{x _ y) 



4(4tt 2 ) 3 J * p 4 

We are left with an IR divergent expression, which is renormalized using the identity 
( |2.4p (with am = 0). Integrating by parts and Fourier transforming back into x space we 
finally find 

I°r = ^^y 2 [ln 2 (x - yfM 2 IR + 2\n{x - y) 2 M 2 R (l - ln(x - y) 2 M 2 ) + b IR ] (4.11) 
Using this result we readily obtain, with k = l/(47r) 4 , 

< 6 a„(-/tA(x-!/)ln(i-j/) 2 M 2 + A(i-j/)/°(i-j/)) , (4.12) 

which, through the use of Bianchi identities V a W a = V^Wa, can be written as an F 
term: 



D x + D 2 = - I d 4 xd 4 yd 2 6 W a (x, 9)W a (y, 9) x 

u{kA(x - y)\n(x-y) 2 M 2 - A(x - y) I° R (x - y)) . (4.13) 



10 



To calculate D$ + D$ it is convenient to decompose the last term of D% (which has an 
overall divergence) into trace and traceless parts: 

4 [(d^A)I°] R = [((4«9 M <% - S fll/ n)A)I°] R + V [(DA)I°] R + c6^6(x - y) . (4.14) 

According to CDR, we have included a local term that can appear in the trace-traceless 
decomposition at the renormalized level [23|. We shall fix its coefficient later on requiring 
gauge invariance. Up to this term, the trace cancels the complete diagram D^. The rest 
of D% is overall UV finite. Renormalizing the sub divergences we find, after some algebra, 



D 3 + D 5 



d*xd*yd*6 T^(x,9)T u (y,9) 



(d„d v - <Vn) - TT AO - y) 1° (a; - y) + k A(x - y) 



ln(x - y) 2 M 2 + 



+ 6. 



K 
U 



□ (A(x - y) \n(x - y) 2 M 2 ) +ed(x- y) 



2 

(4.15) 



With e 



|/c + §• The non-trasverse pieces in the last line must be cancelled by D^. Let 
us consider this contribution next. Again, it is convenient to split it into traceless and 
trace parts: 

T 2 1 T l u V t1\ ( T 2 t2\ , "M" / rl j2 



7-1 



- 6 ~fn - (£, - 5 -fl 2 ) + - ll PP ) + c\J(x - y) 



(4.16) 

We have included again an arbitrary local term in the trace-traceless decomposition. The 
trace can be computed using Gegenbauer polynomials [|13|, [14]]. (Alternatively, its scale 
dependent part can be easily obtained with "systematic" DiffR p7| .) The traceless part 
is UV and IR finite. Therefore, it is fixed by dimensionality and by the traceless condition 
to be of the form a{d,j,d v — ^pd) ■ The coefficient a may be determined from a rather 
cumbersome calculation with Feynman parameters. Adding trace and traceless parts we 
find 



1 Ftfiu 



l Rfiu 



a(d„d„ - 8^D)A(x - y) (4.17) 
- 6, u n A(x - y) In (x - y) 2 M 2 - (|c(3)« + d + ^a)j 



Adding this to (|4.15|) we obtain 
D 3 + D 4 + D 5 = -i J d 4 xd 4 yd 4 9 r»(x, 6)T v {y, 9) 



x (dpd,, - 5^U) [A(x - y) - y) - 2 /dn(x-j/) 2 M 2 - a) 

+ d 4 xd 4 yd A 9 T^x,9)T u (y,9)5 IMU nA , (4.18 
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with h = + C(3)) + fa + c + d . The M- and M/ij-dependent parts are automatically 
transverse, and the non transverse contribution of the scale independent local part can be 
set to zero (so that the sum is gauge invariant) adjusting c and d appropiately. Thanks 
to its transversality, we can then rewrite this expression in a form proportional to the 
classical action: 

D 3 + D A + D 5 = ~ [ d 4 xd A yd 2 9 W a {x,6)W a {y,6) (4.19) 



4 . 
x □ 



A(x - y) (l%,(x - y)) - 2«ln(x - y) 2 M 2 - a) 



Summing up all diagrams and using ( 4.11 ), we finally obtain for the full 2- loop contribu- 
tion in the Feynman gauge (£ = 0) 



E D * = / d 4 xd*yd 2 9 W a (x, 6)W a (y, 9) □ \a(x - y) (4.20) 

i=l 71 J 

x (\n 2 (x-y) 2 M 2 IR + 2\xi(x-y) 2 M 2 IR (l -\n(x-y) 2 M 2 ) + 41n(x-y) 2 M 2 ) 

A possible local and scale-independent contribution has been cancelled by an adequate 
choice of bj R in ( [4.1 1[ ). We see that only single logarithms of M appear in the final 
result, as required by renormalization group invariance (see Section |5|). This fact can also 
be understood in the following way: For a double M logarithm we must have both UV 
sub divergences and overall UV divergences. By power counting, only the TT terms may 
contain overall UV divergences. Furthermore, the traceless parts multiplying IT are finite 
and can have only single logarithms of M, arising from the sub divergences. But gauge 
invariance, i. e. transversality, forces the trace part to have the same logarithm structure, 
so double logarithms of M are also forbidden in the complete result. The same argument 
shows that were the theory scale independent (i.e., finite) to n loops, the background 
two-point function would also be scale independent to n + 1 loopsf]. This line of reasoning 
can be pushed even further if we distinguish in our calculation the two-loop UV scale, M', 
from the one-loop scale, M. Requiring transversality on M and M' independently we see 
that M 1 must cancel in the final result. (This means that the coefficient c contains the 
expression In M' 2 /M 2 .) So, it is only the one-loop scale that appears in the renormalized 
two-loop function. The significance of this observation is discussed in Section |[ 



D The corresponding properties for the 1/e poles in DRed were found in ||. 
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5 Renormalization group equation 



After adding up the partial results (4.7) and ( 4.20Q , the final renormalized expression for 
the background two-point function reads 



G (2)( 



X) 



3C A \n(x 2 M 2 ) 

+ -T^ D 7, " 



+ 



4 2 (47T 2 ^ 2 

'A 



X, 



A 2 (An 



2\2 



X, 



3g 2 C 2 



□ 



ln 2 (a; 2 M 2 R ) + 21n(x 2 M / 2 K )(l - ln(x 2 M 2 )) + 41n(x 2 M 2 ^ 



4 4 (4tt 2 ) 3 
+ 0(0 \ + 0(g 4 ) 



x- 



(5.1) 



Due to background gauge invariance, W undergoes no wave function renormalization. 
Thus the renormalization group equation for this Green function reads 



X , 



. 



(5.2) 



Note that we can only go into the Feynman gauge after evaluating the derivative with 
respect to £. We have not included a term jiRMmd/dMm because in deriving the 
renormalization formula ( |2.4|) we required that the UV and IR scales were independent 
from each other, and hence = M/MiRdMm/dM = 0. In fact, Mm parametrizes 
nonlocal contributions, which are not object of UV renormalization. 

We solve the renormalization group equation perturbatively to order g 2 (two loops). 
The first coefficient in the expansion of 7§(<?) can be read off from the 1-loop vacuum 
polarization for the gauge superfield field V, which is calculated in the Appendix with the 
result 7^(p) = — ^^ 9 2 + 0(g 4 ). With this input, all the nonlocal and scale dependent 
pieces cancel out in (|5.2|) . This is a check of the consistency of our renormalization 
procedure. In particular, note that the IR scale generated in the gauge parameter at two 
loops is cancelled by corrections to the gauge parameter at one loop. The remaining local 



parts of the renormalization group equation (|5.2|) uniquely fix the first two coefficients of 
the beta function: f\ 

2 



0(9) 



3 Ca 3 
'48tt 25 



C a 



9 5 + 0(g 7 ) 



(5.3) 



6 A word on normalization. The coupling constant g as given in ([O]) is V2 times larger than the 
standard Yang-Mills coupling, gym (see p. 55 in Q). In terms of the latter, P(gYAi) 

3/Ca 



3 Ca „2 
' 2 8tt 2 SYM 



3 ( £a\2 5 
'.y^J 9YM 



which matches the expansion of expression ( 1.1 
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6 Discussion 



Using DiffR, we have calculated the complete renormalized correlation function of two 
background gauge superfields in pure N=l SYM. This calculation illustrates the power and 
simplicity of this method in applications to supersymmetric gauge theories. In particular, 
we have seen that DiffR can be employed to subtract IR divergencies as well, and that 
the corresponding IR scale can be clearly distinguished from the UV one. From the 
dependence of the two-point function on the renormalization scale we have derived the first 
and second coefficients of the beta function, b and b\. Furthermore, we have presented 
a new argument showing that the n + 1-loop coefficient vanishes for any supersymmetric 
theory which is finite to n loops. This important property had been proven before using 
DRed §. 

It is interesting to have a closer look at the way in which the two-loop coefficient of 
the beta function, & 1; is generated in our calculation. 

1. UV one-loop sub divergencies are subtracted. This entails a one- loop wave-function 
renormalization of the quantum gauge superfield. The corresponding renormalized 
subdiagrams depend on the one-loop renormalization scale M in a local way. 

2. The overall UV divergencies are subtracted and a new renormalization scale M' 
is introduced. However, the combination of supersymmetry, gauge invariance and 
power counting implies that M' cancels out in the complete renormalized two-point 
function. On the other hand, there remains a nonlocal dependence on M (see Eqs. 
(|1D and flyj))|] 



3. After integration over half the supercoordinates, the dependence on M becomes 
local. 

4. Finally, this local scale dependence is compensated by b\ in the renormalization 
group equation. The off-shell IR divergencies only play a passive role, as they 
exactly cancel in the renormalization group equation. 

Summarizing, the scale associated to the one-loop renormalization of the quantum super- 
field is the one that gives rise to the two-loop coefficient of the beta function! This is 
somewhat surprising because naively one would think that the two-loop coefficient should 
have its origin in M', which is the scale associated to two- loop superficial divergencies. 



7 Note that the expression "dependence on M" refers to the derivative with respect to M and not to 
the term in which M appears, which is always nonlocal. 
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Should this be the case, the beta function would be purely one loop (remember that no 
overall scale can appear at any order n > 1). However, we have seen explicitly in a two- 
loop calculation that the sub divergences play a nontrivial role. More generally, we expect 
that subdivergencies are responsible for all higher-order coefficients of the beta function. 
This agrees with the NSVZ form of the beta function. The fact that M' disappears in 
our method is directly related to the observation in || that in invariant four- dimensional 
regularization methods there are no divergencies after all subdivergencies have been sub- 
tracted. As we have seen, this does not imply b\ = 0. Therefore, it seems that the 
naive perturbative derivation of the beta function from renormalization constants needs 
some modification in this case. This modification is surely related to the presence of the 
anomaly discussed in || || |28|, p9j] . From this point of view, the fact that the standard 
derivation from renormalization constants works in DRed seems related to the fact that 



there are no rescaling anomalies in this method [22 



As a matter of fact, the mechanism we have just described agrees with previous calcu- 
lations in which the corrections to the one-loop result arise from a one-loop anomaly ||, 



U |29|. This anomaly manifests itself in different ways: as a nonzero expectation value 
for the operator WW |§; as the noninvariance of the measure under the rescaling of the 
gauge field [[J; or as the quantum breaking of either supersymmetry or holomorphy in 
the framework of local coupling [B^, |29| . In our description, the anomaly is to be asso- 
ciated with the external loop, and is responsible for the promotion of the M dependence 
into a non- vanishing nonlocal structure that eventually generates b\. This is completely 
analogous to the explicit calculations in SQED performed in || (except for the fact that 
we subtract the subdivergencies). Note that even though M' cancels out, the presence 
of a oo — oo UV behaviour is crucial for the anomaly to exist. On the other hand, as 
emphasized in || , this nonlocal structure is nonanalytic at vanishing external momentum. 
Such on-shell IR divergence, which cancels after integration over the 9 coordinates, is a 



manifestation of the "IR side" of the anomaly []3D| and should not be confused with the 
off-shell IR divergences that we have renormalized in our calculation. More generally, IR 
effects are known to be responsible for quantum corrections to F terms in the 1PI effective 
action [3I|. The on-shell IR divergence arises from the region of virtual momenta of order 
the momentum of the external field. If this region is excluded, one finds b\ = 0. Our 
subtraction of IR divergencies, on the other hand, is local in momentum space and does 
not modify the analiticity properties of the two-point function. The same applies to the 
IR subtractions in ||. Summarizing, the two-loop coefficient of the beta function arises 
from a one-loop UV scale which survives at two loops only when IR effects are included. 
For completeness, we discuss in the rest of this section the connection between the 
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Gell-Mann-Low beta function we have computed (1PI beta function) and the flow of the 
coupling in the Wilsonian action (Wilsonian beta function). The Wilsonian effective ac- 
tion can be understood as the generating functional of Green functions with an IR cutoff, 
for external momenta smaller that this cutoff p3| . Therefore, at least in perturbation 
theory, the ( "holomorphic" ) Wilsonian coupling obeys a one-loop renormalization group 
flow. (For pure Super Yang-Mills this result also holds when nonperturbative effects are 
included.) The higher-order contributions to the physical beta function appear when cal- 
culating expectation values of the operators in the Wilsonian action ||. On the other 
hand, in [g] the "canonical" Wilsonian coupling was shown to obey instead an NSVZ 
flow. According to this work, in a Wilsonian setup this comes about because a rescaling 
of the gauge superfields at each scale is needed in order to keep kinetic terms canonically 
normalized. This rescaling induces an anomalous Jacobian, and it is this anomaly that 
induces the corrections to the one-loop result 0. In this reference, the anomaly is calcu- 
lated a la Fujikawa. A basic element of the calculation is the introduction of a UV cutoff 
and one might believe that the anomaly (and thus the running of the canonical coupling) 
depends only on the UV properties of the theory. However, a closer look shows that the 
low-energy degrees of freedom play a fundamental role |I| [3(J. In fact, the IR degrees 
of freedom must be included in the derivation of the anomaly if low-energy physics is to 
remain unchanged under the field rescalings. In this sense, taking the rescaling anomaly 
into account is equivalent to calculating the expectation value of the Wilsonian action || . 
It was also argued in |J (see |51J as well) that the canonical Wilsonian coupling is 



closely related to the 1PI running coupling. (This implies that the expectation value of the 
canonical Wilsonian action is basically trivial.) This relation has been made more precise, 
in a general context, in f35|. There it is shown that, when all kinetic terms are canonically 
normalized, the Wilsonian coupling becomes independent of the renormalization scale for 
large cutoff. Then one can derive the equation (see also fl36| ) 



Ad/dAg cW (g,A/M) 
[ d/dgg cW (g,A/M) J 



(6.1) 



Here g is the physical coupling, g c w is the canonical Wilsonian coupling, A is the flowing 
cutoff in the Wilsonian action and M is the renormalization scale, introduced by low- 
energy normalization conditions. Using this equation, we show now that (at least for large 
A) the 1PI and the canonical Wilsonian beta functions agree to two loops in perturbation 
theory, as functions of g and g c w, respectively. The flow of g c w is of the generic form 

(3cw(g c w) = Adg cW /dA = b g 3 cW + hg 5 cW + b 2 g 7 cW + 0(g 9 cW ) , (6.2) 
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with constant coefficients b n . On the other hand, g c w can be perturbatively expanded in 
powers of g: 

g cW = g(l + C (A/M)g 2 + C 1 (A/M)g 4 + 0(g e )) , (6.3) 
where we have taken g c w = g at tree level. Inserting ( |6.3| ) into ( |6.2|) we see that 
C (A/M) = 6 l n A + C0; 

C a (A/M) = |6gln 2 A + (26 b cto + 6 1 )lnA + Pl| ( 6 .4) 

where c n are scheme dependent constants. Using (|6.2| ) and (|6.3|) in Eq. (|6.1|) we find 

P(g) = [bo9 3 + hg 5 + (6 2 + 2hC + 3b C 2 - 26 C 1 )( ? 7 + Ofa 9 )]^ , (6.5) 

and using also ( |6.4| ) we see that the beta function is finite: 

P{g) = b g 3 + hg 5 + (b (3c 2 - c x ) + b lCo + b 2 )g 7 + 0(g 9 ) . (6.6) 

Hence, the first two coefficients of the 1PI beta function coincide, in any (mass-independent) 
scheme, with the first two coefficients of the canonical Wilsonian beta function. The other 
coefficients are scheme dependent, as expected. This scheme dependence has been studied 
for general N=l theories in (37]. 
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A Calculation of 7^ 

The classical action in a generic Lorentz gauge reads 

S = -JL J d 8 z (e- 9V D a e 9V )D 2 (e- 9V D Q e 9V ) - (f + 1) J d 8 z (D 2 V)(D 2 V) (A.l) 
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To second order in the quantum field V it reduces to 



S {2) = I fd 8 z VD a D 2 D a V - ^±12 / d s z v(D 2 D 2 + D 2 D 2 )V 



/' 



1 f d 8 z VnU 1/2 v - ^^^ d 8 zVnU v , 



(A.2) 



were we have defined the projectors 



n 



1/2 



n n = 



D a D 2 D a 
□ 

D 2 D 2 + D 2 D 2 
□ 



(A.3) 
(A.4) 



The one-loop correction to this quadratic action can be easily computed. The result in 
the Feynman gauge is 

3C A g 2 



r = 



4(4tt 2 ) 2 
3CV 



d*z l d*z 2 V( Zl )D a D 2 D a V(z 2 )5 12 - - 

{x 1 - % 2 r 

d 8 Zl d 8 z 2 V(z 1 )D a D 2 D a V(z 2 )5 12 D [A(x 12 ) lnx 2 12 M 2 } . (A.5) 



16(4tt 2 

Therefore, the one loop effective action quadratic in the quantum gauge field is 



r = 



d 8 Zl d 8 z 2 V(z 1 )T^ 2 \z 1 ,z 2 )V(z 2 ) 



+ 



d 8 z x d 8 z 2 V(zi) 
3C A g 2 



->)(z 12 )nn 1/2 - { ^^5^(z 12 )m - 



16(4tt 2 



(□[A(^ 2 )lnx 2 2 M 2 ]nn 1/2 + O(0) 



V(z 2 ) . 



(A.6) 



The renormalization group equation for the correlation function of two quantum gauge 
fields has the form: 



To order g 2 (one loop) it is solved for 







(A.7) 



5=0 



Iv 
75 



3CU 2 , 
*9 +■■■ 



8(4tt 

3Ca 2 , 
9 +• 



4(4tt 



(A.8) 
(A.9) 
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